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Abstract 

We study the structure and evolution of neutron stars (NS) the interiors of 
which are modeled using microscopic approaches and constrained by the condi- 
tion that the equation of state (EoS) of matter extrapolated to high densities 
should not contradict known observational data from compact stars and exper- 

■ imental data from heavy-ion collisions (HIC). We use modern cooling simula- 
^ | tions to extract distributions of NS masses required to reproduce those of the 

yet sparse data in the Temperature- Age (TA) plane. By comparing the results 
["T | | with a mass distribution for young, nearby NSs used in population synthesis we 

■ can sharpen the NS cooling constraints. 

1 Introduction 

in 

C*~) . Among compact stars denoted here generically as NSs one can theoretically 

distinguish three main classes according to their composition: hadronic stars, 
quark stars (absolutely bare or with tiny crusts), and hybrid stars (HyS). This 
classification is based on the understanding of the EoS of matter under extreme 
conditions which is at present lively debated. Can astronomy of NSs together 
with experimental data from HICs, where the matter exceeds the nuclear satu- 
ration density allow us to identify the true " name" of a given observed compact 
star? To answer this "simple" question one needs to do a very complex job. 
We have to collect experimentally justified models of nuclear matter, extrapo- 
late their properties to sufficiently high densities, make a compact object model 



CZ2 . 

, based on that and return to observations for comparison of the results. Probably 

this cycle is only one loop of iteration like in any perturbation approach. 

To be sure that one moves in the correct direction one needs to define con- 
straints on our modeling from the investigation of high-density behavior of nu- 
clear matter (NM) activated in such an experiments as planned when the new 
accelerator facility (FAIR) at GSI Darmstadt will be constructed. On the other 
hand the new observations in astronomy provide new limits for the mass and the 
mass-radius relationship, surface temperature etc. of compact stars. These data 
constitute stringent constraints on the equation of state of strongly interacting 
matter at high densities, see [I] and references therein. Of course experiments 
such as CBM (=compressed baryon matter) at FAIR or NICA at JINR which 
are dedicated for the investigation of the phase transition from hadronic matter 
to the quark-gluon plasma (QGP) in HIC have a natural interest for the theory 
of NSs and particularly for the question whether quark matter exists in the core 
of a NS. 



1 Lecture notes for the Helmholtz International Summer School on "Dense Matter 
in Heavy-Ion Collisions and Astrophysics", JINR, Dubna, August 21 - September 1, 
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Here in these lectures, we demonstrate that the present-day knowledge of 
hydrodynamical properties of dense matter allows to construct hybrid EsoS with 
a critical density of the deconfincmcnt phase transition low enough to allow for 
extended quark matter cores and stiff enough to comply with the new mass 
measurements of compact stars [2]. As has been argued by Alford et al. [3] 
[3] there are several modern QCD-motivated quark matter EsoS which could 
provide enough stiffness of high-density matter to be not in conflict with the 
new mass and mass-radius constraints. When compared with hadronic EsoS 
it turns out that quark and hadron matter have rather similar hydrodynami- 
cal properties in the region of the deconfinement transition. This means that 
hybrid stars can masquerade as neutron stars once the parameters of a generic 
phenomcnological quark matter EoS have been chosen appropriately [H [5] • 

However, we suggest a new tool for "unmasking" the composition of neutron 
star interiors which is based on the fact that the state of matter at high densities 
determines the statistics of both NS observables, the temperature-age (TA) data 
as well as the mass distribution. 



2 Equation of state of neutron star matter 
2.1 Hadronic EoS 

There are numerous comparative studies of NM approaches for HIC and NS 
physics applications in which a representation of the NM EoS has been employed 
which is based on the nucleonic part of the binding energy per particle given in 
the form 

E(n,l3)=E (n) + /3 2 E s (n), (1) 

where [3 — 1 — 2x is the asymmetry parameter depending on the proton fraction 
x = rip/n with the total baryon density n = n n + n p . In Eq. (JT|) the func- 
tion Eo(n) is the binding energy in symmetric NM (SNM), and Es{n) is the 
(a)symmetry energy, i.e. the energy difference between pure neutron matter and 
SNM. Both contributions Eo(n) and Es(n) are easily extracted from a given 
EoS for the cases f3 = and (3 — 1, respectively. The parabolic interpolation 
has been widely used in the literature, see e.g. [5], and appears rather robust. 
From Eq. (TTJ all zero temperature EsoS of nonstrange NM can be derived by 
applying simple thermodynamic identities [BJ. In particular, we obtain 

e B (n,(3) = nE{n 7 f3), (2) 

P B {n,fi) = n 2 -^-E(n,/3), (3) 

d \ „ , , / „n _ „9<9 



= \l + n—\E {n)-\3 z TW-p z n—\Es{n) (4) 

for the baryonic energy density s(n) and pressure P(n) as well as the chemical 
potentials of neutron /i n (upper sign) and proton \x p (lower sign), respectively. 

The neutron star matter (NSM) has to fulfill the two essential conditions of 
/3-equilibrium 

A*n = M P + Me = Mp + Mm . ( 5 ) 
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and charge neutrality 



rtp - n e - rip = , 



(0) 



where /_t e and /i^ are the electron and muon chemical potentials, conjugate to 
the corresponding densities n e and n^. In case of a first order phase transition 
a mixed phase could arise in some density interval, see [7] . In general, the local 
charge neutrality condition could be replaced by the global one. However, due 
to the charge screening for the phase transition to quark matter this density 
interval is essentially narrowed [5] • The effect of the mixed phase on the EoS 
is also minor. 

Due to Eq. (jJJJ) the chemical potentials for muons and electrons are equal, 
jUju = /i e . The threshold of muons to appear in the system is determined by 
their mass. In NSM the baryonic and the leptonic parts are considered as an 
ideal mixture in /3-equilibrium 



Under NS conditions one parameter such as the baryo-chemical potential /i^ 
or the baryon density conjugate to it is sufficient for a complete description. In 
/3-equilibrated neutron matter the chemical potential of the neutrons equals the 
baronic one, /i„ = fib- Applying Eq. (J4j) and Eq. (O shows that the electron 
and muon chemical potential can be written as an explicit function of baryon 
density and asymmetry parameter, 



Both electrons and muons are described as a massive, relativistic ideal Fermi 
gas. 

With the above relations only one degree of freedom, namely the baryon 
density, remains in charge neutral and /3-equilibrated NSM at zero temperature. 
Within this description actual properties of NM depend on the behavior of E (n) 
and Es(n) only. Both can be deduced easily from any EoS introduced in these 
lectures. 

Here in the discussions we use the nuclear EsoS that originate from relativis- 
tic descriptions of NM, because in the range of densities relevant NSs and HICs 
relativistic effects are important. 

A brief description of different approaches is given in the following. For more 
details, seepQ and references therein. 

• Phenomenological models 

originate from [T^ and are based on a relativistic mean-field (RMF) de- 
scription of NM with nucleons and mesons as degrees of freedom. 

The mesons couple minimally to the nucleons. The coupling strengths are 
adjusted to properties of NM or atomic nuclei. A scalar meson (a) and 
a vector meson (lu) are treated as classical fields generating scalar and 
vector interactions. The isovector contribution is generally represented by 



s(n,P) 
P(n,0) 



e B (n, (3) + s e (n, (3) + e^(n, /3) , 
P B (n,0)+P e {n,0)+P M (n,l3) . 



(7) 
(8) 



fi e (n,(3)=i(3E s (n). 



(9) 
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a vector meson p. In order to improve the description of experimental data, 
a medium dependence of the effective interaction has to be incorporated 
into the model. 

In many applications of the RMF model, non-linear (NL) self-interactions 
of the a meson were introduced with considerable success and were later 
extended to other meson fields (see ill). As an alternative to NL RMF 
models, appoaches with density-dependent nucleon-meson couplings were 
developed (see [Hj and references within) , where for a flexible description 
of the medium dependence several parameterizations were introduced. In 
the NLp the isovector part of the interaction is described, as usual, only 
by a p meson. The set NLp<5 also includes a scalar isovector meson 5. 

The density dependent RMF models are also represented here by [T3] 
(DD models). They are obtained from a fit to properties of finite nuclei 
(binding energies, radii, surface thickencsses, neutron skins and spin-orbit 
splittings). In the D 3 C model additional couplings of the isoscalar mesons 
to derivatives of the nucleon field are introduced. 

• Dirac-Brueckner-Hartree-Fock (DBHF) approach 

One of the microscopic approaches starting ah initio from a given free 
nucleon-nucleon interaction that is fitted to nucleon-nucleon scattering 
data and deuteron properties is the DBHF method [14] • In this approach 
the nucleon inside the medium is dressed by the self-energy £ based on 
a T-matrix. The in-medium T-matrix which is obtained from the Bethe- 
Salpeter equation plays the role of an effective two-body interaction which 
contains all short-range and many-body correlations in the ladder approx- 
imation. It is possible to extract the nucleon self-energies from DBHF 
calculations which can be compared with the corresponding quantities in 
phcnomenological RMF models, but this is not completely unambiguous 
as discussed in Ref. [15] . Here, we use recent results of (asymmetric) NM 
calculations in the DBHF approach with the rclativistic Bonn A potential 
in the subtracted T-matrix representation. 

• A bridging approach between fully microscopic and more phenomenological 
descriptions 

that adjusts the parameters of the RMF model to results extracted from 
microscopic approaches. As an example of this method, we use a non- 
linear RMF model (KVR) (or slightly modified parameter set (KVOR)) 
with couplings and meson masses depending on the a- meson field 16 . 
The parameters were adjusted to describe the SNM and NSM EoS of the 
Urbana-Argonne group [17] at densities below four times the saturation 
density. In these models not only the nucleon but also the meson masses 
decrease with increasing NM density. Being motivated by the Brown- Rho 
scaling assumption, see |18j . and the equivalence theorem between differ- 
ent RMF schemes, these models use only one extra parameter compared 
to the standard NL RMF model (NL model). 

The sets of parameters and references of used models are given in [T] . 
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Figure 1: The energy per nucleon in SNM Eo(n) (left panel), the symmetry 
energy Eg (n) (middle panel) and the energy per nucleon in NSM (/3-equilibrated 
and charge neutral) for the investigated models (right panel). 

The variation in the NM parameters is directly reflected in the behavior of 
the energy per nucleon in SNM Eo(n) and of the symmetry energy Eg(n) at 
densities above saturation as shown in Fig. [T] The various models of this study 
predict considerably different values for Eq(ji) and Es(n) at high densities. 
Under the condition of /^-equilibrium, however, the range of binding energy per 
nucleon E(n,f3) shows a much smaller variation than expected from Eo(n) and 
Es(n). This is shown in the right panel of Fig. [T] 

2.2 Quark matter EoS 

Since there are no solutions of QCD for quark matter in the nonperturbative 
domain close to the chiral/ deconfinement transition at zero temperature and 
finite density one needs an effective model for the quark matter EoS satisfy- 
ing at least the symmetry requirements of QCD. In order to discuss whether 
dcconfmed quarks can exist in neutron stars a successful class of approaches is 
that of Nambu-Jona-Lasinio (NJL) type models having a Lagrangian with chi- 
ral symmetry which is dynamically broken in the nonperturbative vacuum. In 
our applications in astrophysics we use this microscopic approach, see 1 19) and 
references therein. 

In present investigation we took a three-flavor chiral quark model with self- 
consistently determined quark masses and pairing gaps similar to the paral- 
lel developments in Refs. [211 [22], but generalizing [20] by including an isoscalar 
vector meson current [2] similar to the Walecka model for nuclear matter. While 
a stronger scalar diquark coupling leads to a lowering of the phase transition 
density, an increase in the vector coupling makes the quark matter EoS stiffer. 
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The thermodynamics of the deconfined quark matter phase as described 
within a three-flavor quark model of Nambu-Jona-Lasinio (NJL) type can be 
investigated systematically by applying techniques of finite-temperature field 
theory }23] . The path-integral representation of the partition function is given 

by 



Z(T,A) 



VqDq exp 




q(i@- 



rh + fvy ) q + 



(10) 



= G, 



_a=0,3,8 



-r/D ^2 (qi~f 5 T A \ A Cq T )(q T iCj 5 T A \ A q), 



A=2,5,7 



(11) 



where jl and rh — diagj(m„, ma, m s ) are the diagonal chemical potential and 
current quark mass matrices. For a = 0, tq — (2/3) 1 / 2 ]./, otherwise r a and 
A a are Gell-Mann matrices acting in flavor and color spaces, respectively. C — 
ij 2 j° is the charge conjugation operator and q = <^7°. Gs, rjv, and tjd deter- 
mine the coupling strengths of the interactions. 

The Lorentz three-current, q^q, vanishes in the static ground state of matter 
and is therefore omitted in the vector interaction term in (|10[) , for further details 
see [3|20]. 

After bosonization using Hubbard-Stratonovich transformations, we obtain 
an exact transformation of the original partition function (|10p . In mean-field 
approximation, when the bosonic functional integrals are omitted and the col- 
lective fields are fixed at the extremum of the action we obtain 

Q.MF (T, pi) = - -r- In Z M f (T, fi) 



(3V 
1 



8G.c 



V {m*~m t ) 2 ~— (2wg + $) + — \ Aaa \ 
Vv Vd 

i=u,d,s 1 1 ,4=2,5,7 



-E a /T 



ill — ^0 



(12) 



Here, f2; is the thermodynamic potential for electrons and muons, and f2o is a 
divergent term that is subtracted in order to get zero pressure and energy density 
in vacuum (T = fi = 0). E a (p) are the quasiparticle dispersion relations. Here, 
A AA are the diquark gaps, rh* is the diagonal renormalized mass matrix and 
jl* the renormalized chemical potential matrix, pb* — diagj(/i u — GsrjvWo, [id — 
GsVvUq, /j, s — Gsr)v4>o)- The gaps and the renormalized masses are determined 
by minimization of the mean-field thermodynamic potential (|12[) , subject to 
charge neutrality constraints which depend on the application we consider. In 
the (approximately) isospin symmetric situation of a heavy-ion collision, the 
color charges are neutralized, while the electric charge in general is non-zero. 
For matter in /3-equilibrium, also the electric charge is neutralized. Further 
details are given in Refs. [2] [20 ] [2T ] [22] . 
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Figure 2: Phase diagram for isospin symmetry using the most favorable hybrid 
EoS of the present study. The NM-2SC phase transition is almost a crossover. 
The model DD-F4 is used as a finite-temperature extension of DBHF. For the 
parameter set (rjo = 0.75, r\v = 0.0) the flow constraint is fulfilled but no stable 
hybrid stars are obtained. 



We consider rjo as a free parameter of the quark matter model, to be tuned 
with the present phenomenological constraints on the high-density EoS. Simi- 
larly, the relation between the coupling in the scalar and vector meson channels, 
r]v, is considered as a free parameter of the model. The remaining degrees of 
freedom are fixed according to the NJL model parameterization in Table I of 
[24] . where a fit to low-energy phenomenological results has been made. 

For our investigations of hybrid star structure we use two hadronic matter 
EsoS, where DBHF is an ab-initio calculation for the Bonn-A nucleon-nucleon 
potential within the Dirac-Brueckner-Hartree-Fock approach [T5], discussed in 
the context of compact star constraints in [H [21] . DD-F4 denotes a relativis- 
tic mean-field model of the EoS with density-dependent masses and coupling 
constants adjusted to mimick the behavior of the DBHF approach (13) . The 
transition to the quark matter phase is obtained by a Maxwell construction, 
where the critical chemical potential of the phase transition is obtained from 
the equality of hadronic and quark matter pressures. A discussion of the relia- 
bility of the Maxwell construction for the case of a set of conserved charges is 
given in [8]. The resulting phase diagram for the caseof symmetric NM, relevant 
for applications in HIC experiments as CBM @ FAIR, is shown in Fig. [5J 
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2.3 Constraints on symmetric matter from HIC 



The flow data analysis of dense SNM probed in HICs [25] reveals a correlation 
to the stiffness of the EoS which can be formulated as a constraint to be fulfilled 
within the testing scheme introduced here. 

The flow of matter in HICs is directed both forward and perpendicular 
(transverse) to the beam axis. At high densities spectator nucleons may shield 
the transversal flow into their direction and generate an inhomogeneous density 
and thus a pressure profile in the transversal plane. This effect is commonly 
referred to as elliptic flow and depends directly on the given EoS. An analysis 
of these nucleon flow data, which depends essentially only on the isospin in- 
dependent part of the EoS, was carried out in a generic model in Ref. [55]. In 
particular it was determined for which range of parameters of the EoS the model 
is still compatible with the flow data. The region thus determined is shown in 
Fig. [3] as the dark shaded region. Ref. [H] then asserts that this region limits 
the range of accessible pressure values at a given density. For our purposes we 
extrapolated this region by an upper (UB) and lower (LB) boundary, enclosing 
the light shaded region in Fig. [3] 

Its upper boundary is expected to be stable against temperature variations. 
The important fact is that the flow constraint probes essentially only the sym- 
metric part of the binding energy function Eo(n). 

In the hybrid star two phases are separated with a jump in baryon density 
corresponding to derivative of the pressure with respect to the baryochemical 
potential. The same is shown also for isospin-symmetric matter in Fig. [3] A 
slight variation of the quark matter model parameters r\o and r\y results in con- 
siderable changes of the critical density for the phase transition. The appearent 
problem of a proper choice of these parameters we solve by applying the flow 
constraint shown in FigJ3] At first we fix the vector coupling by demanding 
that the high density behavior of the hybrid EoS should be as stiff as possible 
but still in accordance with the flow constraint. We obtain r\y — 0.50 such 
that the problem of the violation of the flow constraint for the DBHF EoS at 
high densities is resolved by the phase transition to quark matter. The optimal 
choice for t\d is thus between 1.02 and 1.03. In the following we will investi- 
gate the compatibility of the now defined hybrid star equation of state with CS 
constraints. 

3 Compact star structure, Mass- Radius relation 

The mass and structure of spherical, nonrotating stars, to which we limit our- 
selves in this paper, is calculated by solving the Tolman-Oppenheimer-Volkov(TOV)- 
equation, which reads as 



dP(r) 
dr 



G[e(r) + F(r)][m(r) + 47rr 3 P(r)] 
r[r — 2Gm(r)] 



(13) 
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Figure 3: Pressure vs. density region consistent with experimental flow data 
in SNM (dark shaded region) . The light shaded region on the left panel extrap- 
olates this region to higher densities within an upper (UB) and lower border 
(LB). A set of hadronic EoS is shown for comparison in the left panel, see text 
and Ref . [T] . The right panel shows a comparison to EsoS with a deconfinement 
phase transition from DBHF hadronic matter to superconducting quark matter 
with isovector mean field taken from Ref. 2J. The quark matter EoS favored 
by the flow constraint has a vector coupling r\y — 0.50 and a diquark coupling 
between r]D = 1-02 (blue solid line) and r\n = 1.03 (black fat solid line); results 
for four intermediate values r\D — 1.022 . . . 1.028 are also shown (thin solid 
lines). 

where the gravitational mass m(r) inside a sphere of radius r is given by 



which includes the effects of the gravitational binding energy. Eq. (fLUf describes 
the gradient of the pressure P and implicitly the radial distribution of the energy 
density e inside the star. The EoS, i.e., the relation between P and e leaves 
as much unknowns as we need to solve the set of differential equations. We 
supplement our EsoS describing the NSs interior by an EoS for the crust. For 
that we use a simple BPS model [57] • Due to uncertainties with different crust 
models one may obtain slightly different mass-radius relations. 

The stellar radius R is defined by zero pressure at the stellar surface, P(R) = 
0. The star's cumulative gravitational mass is then given by M = m(R). Each 
configuration is specified by the given central density value n(0) at radius r = 0. 



r 




(14) 



o 
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Figure 4: Left panel: Mass-Radius constraints from thermal radiation of the iso- 
lated NS RX J1856.5-3754 (grey hatched region) and from QPOs in the LMXBs 
4U 0614+09 (green hatched area) and 4U 1636-536 (orange hatched region) 
which shall be regarded as separate conditions to the EsoS. For the NS in 4U 
1636-536 a mass of 2.0 ± 0.1 Mq is obtained, so that the weak QPO constraint 
would exclude the NLp and NLpS EsoS whereas the strong one would leave only 
DBHF, DD and D 3 C. Right panel: Mass - radius relationship for CS sequences 
corresponding to a nuclear matter EoS (DBHF) and different hybrid star EsoS 
(DBHF+NJL), see text and [2]. Indicated are also the constraint on the mass 
from the pulsar J0751+1807 [28 and on the mass-radius relationship from the 
isolated neutron star RX J1856 [53] The dotted lines indicate the gravitational 
redshift, z = (1 — 2GM/R)~ 1 / 2 — 1, of photons emitted from the compact star 
surface. Present constraints on the mass-radius relation of CSs do not rule out 
hybrid stars. 

3.1 Observational constraint on the Mass and Radius 

On the plot of the Mass-Radius relationship in Fig. 0] we also show regions 
corresponding to different observed objects. These results constrain the model 
EsoS for both hadronic and hybrid star matter. These data are the following 

• Mass-Radius relation from LMXBs 

as implies [30] IS] the observation of a source of quasi-periodic brightness 
oscillations (QPOs) with a maximum frequency is f max limits the stellar 
mass and radius to 

M < 2.2 M Q (1000 Hz/zw)(l + 0.75j) , . 

R < 19.5 km(1000 Hz/z/ max )(l + 0.2j) , 1 ' 

where j = cJ/GM 2 ~ 0.1 — 0.2 is the dimensionless stellar angular mo- 
mentum. 
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• Mass- Radius relation from RX J 1856: 

The thermal radiation of nearby isolated NS RX J1856. 5-3754 with black- 
body spectrum temperature = 57 eV [35] because of initially estimated 
distance (60 pc) has been considered as a self-bound strange quark star 
with radius Roq m 8 km. For the new measurements (distance ~ 117 pc), 
RX J1856 has to have a rather large radius of R ~ 14 km (Roo = 16.8 
km), when the mass is 1.4 M© [55] , 

• Maximum mass constraint 

from measurements on PSR J0751+1807 imply a pulsar mass of 2.1 ± 
0.2 (1q 5) (first error estimate with la confidence, second in brackets 
with 2a confidence) [2"5] , 

The resulting lower bound in the mass radius plane is shown in Fig. H] 
Particularly interesting is the constraint from RX J 1856 and there are three 
ways to interpret this result: 

A) RX J1856 belongs to compact stars with typical masses M ~ 1.4M© and 
would thus have to have a radius exceeding 14 km (see Fig. H]) . None of 
the examined EsoS can meet this requirement. 

B) RX J1856 has a typical radius of R ~ 12 — 13 km, implying that the EoS 
has to be rather stiff at high density in order to allow for configurations 
with masses above ~ 2 M©. In the present work this condition would be 
fulfilled for DBHF, DD and D 3 C. This M > 1.7 M© explanation implies 
that the object is very massive and it is not a typical NS which would 
have M < 1.5 M©, as follows from population synthesis models. 

C) RX J1856 is an exotic object with a small mass ~ 0.2 M©, which would 
be possible for all EsoS considered here. No such object has been observed 
yet, but some mechanisms for their formation and properties have been 
discussed in the literature |33j . 

An other constraint on the EoS of stellar matter could be the critical mass 
Mdu corresponding to the threshold of the most effective cooling mechanism, so 
called Direct Urea (DU) processes n — > p + e~ + P e , which is an ultimate cause 
of very fast cooling of stars when its mass is only slightly above the critical 
value Mdu [S3]- For each hadronic EoS the critical mass is shown in Fig. [5] 
It depends on the symmetry energy E$ and occurs when the proton fraction 
exceeds 11% (without muons and about 14% with muons in the medium) pQ. 

To consider the DU as a strong constraint one needs to investigate the cooling 
problem together with the nucleon superfluidity which suppresses the cooling 
rates. Calculations show that indeed the values of the pairing gaps used in the 
literature are not sufficient to resolve the DU as a problem [35] . 

4 Cooling behavior of hybrid stars 

The cooling behavior of compact stars belongs to the most complex phenomena 
in astrophysics. Therefore, the codes for its numerical simulation as developed 
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Figure 5: Mass versus central density for compact star configurations obtained 
by solving the TOV equations (fl"5|) and (fT4"|) for all EsoS introduced in Subsect. 
2.1. Crosses denote the maximum mass configurations, filled dots mark the 
critical mass and central density values where the DU cooling process becomes 
possible. According to the DU constraint, it should not occur in "typical NSs" 
for which masses are expected from population synthesis [36] to lie in the lower 
grey horizontal band. The dark and light grey horizontal bands around 2.1 M Q 
denote the la and 2a confidence levels, respectively, for the mass measurement 
of PSR J0751+1807 [21]. 



by a few groups contain inputs (cooling regulators) of rather different kind, see, 
e.g., [37l EH1 [39l [4Ql [41] . Attempts to develop a Minimal Cooling Paradigm [38] 
by omitting important medium effects on cooling regulators [401 135] unfortu- 
nately result in inconsistencies and suffer therefore from the danger of being not 
reliable. To develop a paradigmatic cooling code as an open standard, however, 
is rather necessary to cross-check the present knowledge of the groups before 
more sophisticated mechanisms like anisotropies due to the magnetic field [37] 
or special processes in the NS crust or at the surface are taken into account. 
Therefore, it is still premature to attempt an identification of the NS interior 
from the cooling behavior. 

In order to circumvent such a model dependence we employ a given cooling 
code developed in Refs. [391 H2] and vary the matter properties such as EoS, 
superconductivity and star crust model such as to fulfill all constraints known 
up to now (mass, mass-radius, TA, brightness, etc.). Moreover, we try to use 
consistent inputs. 

The main neutrino cooling processes in hadronic matter are the direct Urea 
(DU), the medium modified Urea (MMU) and the pair breaking and formation 
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(PBF) whereas in quark matter the main processes are the quark direct Urea 
(QDU), quark modified Urea (QMU), quark bremsstrahlung (QB) and quark 
pair formation and breaking (QPFB) [33]. Also the electron bremsstrahlung 
(EB), and the massive gluon-photon decay (see [33]) are included. 

The 15*o neutron and proton gaps in the hadronic shell are taken according 
to the calculations by [35] corresponding to the thick lines in Fig. 5 of Ref. 
•10) . However, the 3P2 gap is suppressed by a factor 10 compared to the BCS 
model calculation of [45] , consistent with arguments from a renormalization 
group treatment of nuclear pairing [3B] . Without such a suppression of the 3P2 
gap the hadronic cooling scenario would not fulfill the TA constraint, see |35j . 

The possibilities of pion condensation and of other so called exotic processes 
are included in the calculations for purely hadronic stars but do not occur in 
the hybrid ones since the critical density for pion condensation exceeds that 
for deconfinement in our case [31]. While the hadronic DU process occurs in 
the DBHF model EoS for all neutron stars with masses above 1.27 Mq, it 
is not present at all in the DD-F4 model, see the right panel of Fig. [6] We 
account for the specific heat and the heat conductivity of all existing particle 
species contributing with fractions determined by the (3— equilibrium conditions. 
Additionally, in quark matter the massless and massive gluon-photon modes also 
contribute. 

In the 2SC phase only the contributions of quarks forming Cooper pairs 
(say red and green) are suppressed via huge diquark gaps, while those of the 
remaining unpaired blue color lead to a so fast cooling that the hybrid cooling 
scenario becomes unfavorable 47J. Therefore, we assume the existence of a weak 
pairing channel such that in the dispersion relation of hitherto unpaired blue 
quarks a small residual gap can appear. We call this gap Ax and show that 
for a successful description of the cooling scenario Ax has to have a density 
dependence. We have studied the ansatz Ax = Ao exp [—a (/VMc — 1)], where 
\i is the quark chemical potential, /i c = 330 MeV. For the analyses of possible 
models we vary the values of a and Ao, given in the Table 1 of 48] and shown 
in the left panel of Fig. E] 



The physical origin of the X-gap is not yet identified. It could occur, e.g., 
due to quantum fluctuations of color neutral quark sextet complexes [49 . Such 
calculations have not yet been performed within the relativistic chiral quark 
models. The size of the small pairing gaps in possible residual single color/ 
single flavor channels [50] is typically in the interval 10 keV - 1 MeV, see 
discussion in [51] . The specific example of the CSL phase is analyzed in more 
in detail in Refs. El El El EH] . 

4.1 Cooling curves in the TA diagram 

We consider the cooling evolution of young neutron stars with ages t ~ 10 3 — 10 6 
yr which is governed by the emission of neutrinos from the interior for t < 10 5 yr 
and thermal photon emission for t > 10 5 yr. The internal temperature is of the 
order of T ~ 1 keV. This is much smaller than the neutrino opacity temperature 
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n[fm" 3 ] n(r = 0) [fm" 3 ] 



Figure 6: Density dependence of the pairing gaps in nuclear matter together with 
that of the hypothetical X-gap in quark matter (left). Mass-central density relation 
for the two hadronic EoS models DBHF and DD-F4 (right). The dot indicates the 
onset of the DU process. 

Topac ~ 1 MeV as well as critical temperatures for superconductivity in nuclear 
(T c ~ 1 MeV) or quark matter (T c ~ 1 - 100 MeV). Therefore, the neutrinos 
are not trapped and the matter is in a superconducting state. In Fig. [S]we show 
the density dependence of the pairing gaps in nuclear matter |45l 139] together 
with that of the hypothetical X-gap in quark matter [56l [47j [48] . The phase 
transition occurs at the critical density n c = 2.75 no = 0.44 fm~ 3 . 



In Fig. [7j we present TA-diagrams for two different hadronic models and in 
Fig. [5] for two hybrid star cooling models presented in Ref . [48 . In the cooling 
calculations presented the crust model, i.e. the T m — T s relationship between 
be temperatures of the inner crust and the surface [3JI] , has been chosen as such 
to fulfill the TA test: each data point should be explained with cooling curve 
belonging to an admissible configuration. 

The TA data points are taken from |38] . The hatched trapeze- like region rep- 
resents the brightness constraint (BC) [57]. For each model nine cooling curves 
are shown for configurations with mass values corresponding to the binning of 
the population synthesis calculations explained in [48] . 

In [48] for these hybrid cooling scenario the logN-LogS distribution con- 
straint has been considered and it has been suggested to use the marking of TA 
diagram with five grey values in order to encode the likelihood that stars in that 
mass interval can be found in the solar neighborhood, in accordance with the 
population synthesis scenario, see Fig. [SJ The darkest grey value, for example, 
corresponds to the most populated mass interval 1.35 - 1.45 Mq predicted by 
the mass spectrum used in population synthesis. 



14 




Figure 7: Hadronic star cooling curves for DBHF model EoS. Different lines corre- 
spond to compact star mass values indicated in the legend (in units of Mq), data 
points with error bars are taken from Ref. [38] . 

On the TA- diagram we show also the other constraint given by the maximum 
brightness of CSs, as discussed by |57]. It is based on the fact that despite many 
observational efforts one has not observed very hot NSs (logT > 6.3 — 6.4 K) 
with ages of 10 3 - 10 4 5 years. Since it would be very easy to find them - if 
they exist in the galaxy - one has to conclude that at least their fraction is very 
small. Therefore a realistic model should not predict CSs with typical masses at 
temperatures noticeable higher than the observed ones. The region of avoidance 
is the hatched trapezoidal region in Figures [7] and [HI 

Altogether, the hybrid star cooling behavior obtained for our EoS fits all of 
the sketched constraints under the assumption of the existence of a 2SC phase 
with X-gaps. 

5 Unmasking neutron star interiors using cool- 
ing simulations 

Considering the interior of neutron stars as a " laboratory" where matter under 
conditions of extreme densities occurs [5 8) we hope that it is possible to identify 
the composition of neutron star interiors from their observable properties like, 
e.g., masses, radii, rotational and cooling evolution [5^1 150], 

The question for the existence of quark matter inside a compact object has 
often been mismatched with that of the existence of strange stars, essentially 
made up of strange quark matter without a hadronic shell and thus very compact 
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Figure 8: Cooling curves for hybrid star configurations with 2SC+X pairing pattern 
and X-gap model I (left) versus model IV (right). For the gaps see the left panel 
of Fig. 1. The grey value for the shading of the mass bin areas corresponds to the 
probability for that mass bin value in the population synthesis model of Ref. }3(T. 



with radii less than ~ 10 km [5TT Ifj2"] . The measurements of high masses for 
objects like the pulsar PSR J0751+1807 gE] with M = 2.1 ± 0.2 M Q or of the 
mass-radius relation from the RX J 1856-3754 pointing to either large masses 
or large radii of R > 14 km [29] require a stiff equation of state and exclude 
standard models for hyperonic or quark matter interiors as well as mesonic 
condensates, see [S"5] . 

Therefore, all suggested signals which are based on large first order phase 
transition effects and changes in the mechanical roperties such as the timing 
behavior of pulsar spin-down [64] . frequency clustering [65 or population clus- 
tering H7j of accreting NSs should not be applicable. 

Here we discuss a new, sensitive tool for "unmasking" the composition of 
neutron stars which is based on their cooling behavior [68j . As the cooling 
regulators such as neutrino emissivities, heat conductivity and specific heat in 
quark matter might be qualitatively different from those in nuclear matter, 
due to the chiral transition and color superconductivity with some possibly 
sensible density dependence, the TA curves for hybrid stars could be significantly 
different from those of neutron stars. 

5.1 Mass distribution from TA data 

In order to reach the goal of unmasking the neutron star interior we use here 
a new method for the quantitative analysis of the cooling behavior consisting 
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in the extraction of a NS mass distribution from the (yet sparse) TA data [55] 
and its comparison with the (most likely) mass distribution from population 
synthesis models of NS evolution in the galaxy J35] . 

This method can be described as follows. For a given cooling model defined 
by the EoS and the cooling regulators, to each configuration with a gravitational 
mass M corresponding cooling curve T(t; M) can be determined. For each set 
of mass values Mi, i = 0, . . . , Nm which is defining the borders of Nm mass bins 
as, for example, in the population synthesis one can mark a region with borders 
of pair of neighboring cooling curves T(t; Mi) and T(t; Mj_i). This strip in the 
TA plane corresponds to the i th mass bin. As a measure for the number of 
cooling objects to be expected within this mass bin we chose 



where N coo \ denotes the total number of observed coolers used for the analysis 
and Pj(T,t) is the probability density to find the j th object at the point (T,t) 
in the TA plane. In Ref. [68] the simplest ansatz has been made that Pj(T, t) is 
constant in the rectangular region defined by the upper and lower limits of the 
confidence intervals corresponding to the temperature and age measurements, 
(Tji,T JU ) and (tji,tj u ), respectively, 



P iT,t) = [(^-i i „)(T^-T i „)]- 1 e(T-T J -06(T^-T)e(t-^)efe„-t) . (17) 



Note that in the case when the exact age the object is known (e.g., for a historical 
supernova), the time-dependence of Pj(T,t) degenerates to a 8- function and 
the i-integral in (fTB| can be immediately carried out, leaving us with a one- 
dimensional probability measure. 

This method has been applied to the cooling models for hadronic and hybrid 
stars described in the previous section. The results for the extracted mass dis- 
tributions are normalized to 100 objects, defining N(M) = 100 iVi/Q^i=!o -^)> 
and shown in Fig. [5] 

As we see from Fig. the results are very sensitive to the chosen cooling 
model. In the hadronic scenario the onset of the DU cooling mechanism drasti- 
cally narrows the mass distribution around the critical mass for the DU onset, 
see Fig. [7J On the other hand the slow cooling model predicts more massive 
objects than could be justified from the independent population analysis. 

When comparing the density dependence of the pairing gaps, given in the 
left panel of Fig. El with the extracted mass distributions for the corresponding 
hybrid models in the right panel of Fig. [9l the direct relationship between the 
superconductivity and the mass distribution becomes obvious. 

The DU problem as it was previously discussed in the literature [33 [TJ [5j5] 
was based on the intuitive understanding that the mass distribution can not 
be peaked at a critical mass value which accidentally is unique for all observed 
young objects. Our modification of the definition of the DU problem does not 
contradict that suggestion, but rather provides a quantitative measure which 
to some extent rehabilitates the validity of cooling scenarios including the DU 
process. 




(16) 
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On the other hand, the EoS model should obey the mass constraints too. 
Therefore, using the models discussed in this work we demonstrate that the 
most preferable structure of the compact object is likely to be a hybrid star 
with properly defined color superconductivity of the quark matter state in the 
core. 
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Figure 9: NS mass spectra extracted from the distribution of cooling data for both 
hadronic EoS models (left panel) and for hybrid stars with X-gap models I-IV (right 
panel) . For comparison, the mass disctribution of young, nearby NS from the popula- 
tion synthesis of Popov et al. |48] is shown at the bottom of the panels. 



6 Conclusions 

In these lectures we have discussed the recent developments towards a test 
scheme for the EoS of matter under the extreme conditions using observational 
data from compact stars and experimental results from heavy-ion collisions. 

The application of this scheme to specific EsoS offers some interesting in- 
sights which reveal the discriminative power of their combined tests in a broad 
region of densities including the possibility of a deconfinement phase transition. 

In addition to the astrophysical part of the tests the results for the elliptic 
flow [26] restrict the hardness of symmetric matter which is in opposition to the 
constraint from neutron star maximum mass and/or mass-radius measurements. 

One of the main results concludes that no present phenomenological finding 
bears a strong argument against the presence of a QM core inside NSs [H [2] ■ 
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The resulting phase diagram for stellar as well as symmetric matter includ- 
ing the transition from nuclear to quark matter and satisfying the constraints 
discussed here exhibits almost a crossover transition with a negligibly small co- 
existence region and a tiny density jump. Therefore, the differences between 
the mechanical properties of the hadronic and hybrid stars are very tiny and 
the problem of the cooling history of the compact stars becomes central for the 
discrimination between alternatives for the composition of their interiors. 

We have presented mass distributions obtained by analysing cooling calcu- 
lations for NSs and demonstrated that they provide a sensible measure for the 
composition of compact star matter, even if the mechanical properties of the 
compact objects are almost identical. 
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